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Abstract 

In our preceding serial works, we have investigated the generation of higher- 
order atomic dipole squeezing (HOADS) in a high-Q micromaser cavity, dis- 
cussing the effects of dynamic Stark shift, atomic damping, atomic coherence 
and nonlinear one-photon processes and different initial states (for example, 
correlated and uncorrelated states, superposition states, squeezed vacuum). 
In this paper, we continue to study HOADS in a high-Q micromaser cavity, 
but consider that the atom interacts with the optical field via a multi-photon 
transition process and that the initial atom is arbitrarily prepared. For a 
vacuum initial field, we demonstrate that HOADS cannot occur if the atom is 
initially prepared in a chaotic state and that a coherent atomic state generates 
less efficient and stable HOADS than an arbitrary one. It is found that large 
detuning may lead to enhanced and strong HOADS. 
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1. INTRODUCTION 



In 1926, Heisenberg [1] proposed a fundamental, general principle of nature, known as 
Heisenberg uncertainty principle. Although this principle is of no consequence for the macro- 
scopic world, it does play an important role in dealing with problems met in the microscopic 
world (for example, atom, molecule, nanocrystals, solid, and so on). It was Kennard [2] 
who presented us the first example of nonclassical states, called squeezed state, in which the 
quantum fluctuations in a dynamical observable may be reduced below the standard quan- 
tum limit at the expense of increased fluctuations in its canonical conjugated one without 
violating the Heisenberg uncertain principle. Furthermore, Plebahski [3] made an important 
contribution to the theory of squeezed states in 1956. The exciting and constructive experi- 
mental results of the first successful generation and detection of squeezed states was reported 
in the middle of the 1980's [4]. Since then, the squeezed states have been extensively studied, 
both theoretically and experimentally [5-17], due to their potential applications in quantum 
communication, high- resolution laser spectroscopy measurement, gravity wave detection and 
quantum information theory (quantum teleporation, dense coding, cryptography, quantum 
nondemolition measurement and power-recycled interferometer). These extensive studies 
have shown that a lot of nonlinear optical systems could generate squeezed states for the 
field and atom, for example, in two-photon laser, parametric amplifiers, four-wave mixing, 
resonance fluorescence, Rydberg atom maser and cooperative Dicke system. A general rela- 
tionship between field and atomic dipole squeezing has even been established under different 
initial conditions for both the field and atom [5,8,10-12]. 

Recently, the rapid development of techniques for making higher-order correlation mea- 
surements has resulted in an increased interest in generating another kind of nonclassical 
states, known as higher-order squeezed states [18]. In these states, higher-order quantum 
fluctuations in one quadrature of the field or atomic dipole could be reduced without violat- 
ing higher-order uncertainty relations [18]. It has been predicted that higher-order squeezed 
states for the radiation In 1997, one of us, Xie, with his colleagues [19] introduced the 
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concept of higher-order atomic dipole squeezing (HOADS) and applied it to high-Q micro- 
maser cavities. In our serial works [20-25], we have discussed thoroughly the important 
connection of HOADS with the second-order field and atomic dipole squeezing (ADS), the 
effects of a nonlinear one-photon process and dynamical Stark shift in two-photon processes 
on HOADS, and different initial conditions regarding both the atomic and field (for exam- 
ple, uncorrelated and correlated coherent states, superposition states, squeezed vacuum). 
These extensive studies have made an important contribution to the theory of higher-order 
squeezed states and provided us an approach to extracting information efficiently from an 
optical signal by higher-order correlation measurements. However, the initial conditions are 
more or less individualized, and only one- and two-photon transitions have been investi- 
gated. Hence, it is our purpose in this paper to treat this problem in a more general way. 
As a first step, we consider that the initial atom is arbitrarily prepared and that the atom 
interacts with the field via a multi-photon transition process. In this paper, we use the two- 
level multi-photon Jaynes-Cummings model [9,26,27] and utilize the density matrix method 
which can describe very well arbitrary atomic-field states such as a mixed state. Their actual 
squeezing conditions and behaviours were carefully discussed. 

This paper is organized as follows. In section 2, we describe the theoretical model 
of multi-photon interaction between a cavity field and a two-level atom and introduce the 
concept of HOADS. In section 3, assuming that the atom is initially prepared in an arbitrarily 
state, we investigate the generation of HOADS in the multi-photon Jaynes-Cummings model. 
A summary is given in the last section. 

2. THEORY 

Here we consider the interaction of a two-level atom with a single-mode quantized field 
involving the emission or absorption of multi-photons per atomic transition. The Hamilto- 
nian of this system can be written as follows 

H = na)a + uS z + e(aK + a ? )(S_ + S+) (1) 
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which is the multi-photon Jaynes-Cummings model without the rotating-wave approxima- 
tion (RWA) [9,26]. If the RWA is explicitly used, then we have 



H = Qa^a + uS z + e{a^S- + a^S + ) 



(2) 



which is the multi-photon Jaynes-Cummings model within the RWA [27]. S z and S± are 
operators of the atomic pseudospin inversion and transition, respectively, and satisfy the 
commutation [7]: [S+, S-] = 2S Z and [S z , S±] = ±S±. oo is the transition frequency for the 
atom, at and a are the creation and annihilation operators for the photons with the frequency 
Q, which obey the boson operators' commutation relations, [a, at] = 1. e is the coupling 
constant between the atom and the radiation field, and £ is the absorbing or emitting photon 
numbers per atomic transition. Throughout we employ the unit with h = c = 1. 

In order to investigate the squeezing properties of the atomic dipole variables, we follow 
the standard procedure of defining the slowly varying operators [7] 

1 



2 
1 
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S+e-™ + S_e iwt 



(3) 
(4) 



where S x and S y , in fact, correspond to the dispersive and absorptive components of the 
slowly varying atomic dipole [28], respectively. One can easily show that the above operators 
obey the commutation relation, [S x ,S y ] = iS z . Correspondingly, we found the higher-order 
uncertainty relation [19] concretely given by 



(AS x ) p (AS y ) p > \ \[{AS x ) p /\ [AS y ) p / 2 ] 



(5) 



For a two-level atom, we have 



(A^) p = 2- p + £ (C* - C^)2 k - p < S 3 > k (j = x,y), 

fc=2,4,... 



(6) 



P/2-1 



[(A^) p / 2 , (AS y ) p / 2 ] =i<S z > C^ /2 C n P/2 2 2 - p+m+n < S x > m < S y 



m,n 



(7) 
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where P/2 > 1 is an integer. If P/2 is even (odd), then m and n are odd (even) numbers. 
It is convenient that we define the following functions 

F^P) = (AS X ) P - \ \\^S x ) p l\ (AS y ) p ^]\ , (8) 
F 2 (P) = (AS y ) p - \ \\{AS x ) p l\ (AS y ) p l 2 \\ . (9) 

Then, higher-order quantum fluctuations in the component S x (or S y ) of the dipole are 
squeezed if Fi < (or F 2 < 0). This is the general definition of HOADS [19]. 

3. Results 

We denote | n) as the Fock state of the radiation field and | +) and | — ) as the excited and 
ground states of the two-level atom, respectively. Using the standard bare-state procedure 
introduced by Xie et al. [29], we can calculate the density operator of the atom- field coupling 
system at time t with an arbitrary initial condition p(t = 0) by 

p{t) = u(ty P (t = o)u(t) = ^PjM^ i{Ei ~ Ej)t I Mi I (io) 

ij 

where Pji(0) = (j \ p(0) | i). \ i) and Ei are the ith eigenstate and its corresponding 
eigenvalue, respectively, and U(t) = exp(-iHt) is the unitary evolution operator. Therefore, 
the expectation value of any physical operator O at time t can be arrived at through 

(0(t)) = Tr[p(t)O(0)} = | p(t)O(0) | i) 

i 

oo oo 2 2 

= J2pij(t)U I o(o) | <> = E E EE^)(*^ I o(o) i * n , Q ), (ii) 

ij n=0 m=0 a=l [3 

where | i) = E™=o ELi I ^n,a)- Here | ^ n ,a) denotes the dressed state with n being the 
photon number of the optical field. The wave function | ^ n ,a=i,2) can be expressed by 



^n,i) = (|+,n)+|-,n + p))/ v / 2, (12) 
y n ,2) = (\+,n)-\-,n + p))/V2 (13) 



with the corresponding eigenenergies 



E n>1 = (n + £/2)u + eIT n , (14) 
E nt2 = (n + Z/2)u-eIL n , (15) 



where Il n is defined by 



As usual, we assume that the initial density matrix p(t = 0) can be decomposed into its 
atomic and field parts, i.e., p(0) = p a (0) <8> p/(0). In this paper, we assume that the initial 
atomic state is arbitrarily prepared, i.e., 

p a (0) = sin 2 9 | +)(+ | + cos 2 9 I -)(- I +77 (e"** | -) (+ | +e l<t> | +> (- |) (17) 

with 

t/ 2 < sin 2 # cos 2 # < j, (18) 
and the field is in its vacuum state, i.e., 

P/(0)=|0)(0|. (19) 
Then, we have arrived at the expectation values of S x , S y and S z 

< S x (t) > = 77 jcos(rY) cos + </>)- A^^ sin (|r + «/>) | , (20) 

< 5,(0 > = 77 jcos(rr) sin (|r + «/>)- cos (|r + 0)1, (21) 



1 f • 2 fl A 2 /4 + £! cos(2]>) \ 
<S z (t)> =- sm 2 — i ^ — ^ i - cos 2 I , 



(22) 



where r = et is the scaled interaction time, A = (£fi — uo)/e is the scaled detuning, T 



^A 2 /4 + £! is the scaled Rabi frequency and is the relative phase between the excited 
state I +) and the ground state | — ) of the two-level atom. If we take 77 = sin 9 cos 9, the 
above results are in agreement with those of Xie et al. [9,26]. 

If the two-level atom is initially prepared in a chaotic state (i.e., 77 = 0), one can easily 
see that the second- and higher-order quantum fluctuations in S x or S y cannot be squeezed. 



In the case of pure initial atomic states, which still keep all the relevant physical features, 
we have rj —\ sin (20) | /2 and find that no HO ADS occurs for a completely excited or ground 
atom, but HOADS could appear for a coherent atom as shown below. 

One interesting result is the so-called population trapping [30] , in which the expectation 
value of S z takes on a steady value, i.e., d < S z (t) > / dt = 0. This results in \/£\t = kix or 
(k + l/2)n (k=0, 1, 2, ...). Obviously, for y^Tr = (k + l/2)ir, the second- and higher-order 
quantum fluctuations in S x or S y cannot be squeezed. However, for y^Tr = kn, HOADS 
occurs for proper choice of the parameters rj, 9 and 0. Details are given below. 

In Fig.l, we show a 3D contour plot for the function F\(P — 6), 6 and 77 for = (fc+l/2)7r 
(k=0,l,2,...) in the case of resonance. We find that for < rj < 0.5 the sixth-order quantum 
fluctuations in S x can be squeezed by properly choosing the parameter 9. Especially, when 
77 = 0.5, HOADS appears almost in the whole range of 9 except a small range around 9 = 7r/4 
or 9 = 3n/4. Also it is obvious that the pure atomic state with rj —\ sin (29) | /2 generates 
less HOADS than the atomic state with < rj < 0.5. However, as we vary the phase to 
= kir (k=0,l,2,...), as shown in Fig.2, there exists only a small region of HOADS in the rj-9 
parameter space. This implies that HOADS is much sensitive to the relative phase between 
the excited and ground states of the two-level atom. 

In Fig. 3 and 4, we show 3D contour plots for the function F\(P = 6), parameter 9 
and phase for 77 =| sin(20) | /2 and r/ = 0.5, respectively, for a resonance transition 
(i.e., A = 0). HOADS patterns are clearly shown in the 9-<p parameter space. We find 
that HOADS is much sensitive to the phase for an initial pure atomic state, whereas it is 
relatively stable with respect to the phase for an arbitrary state with 77 = 0.5. This implies 
that an arbitrary atomic state could generate much efficient and stable HOADS than a pure 
one. 

In Fig. 5, we show the time evolution of the function -F\(P = 6) for 9 = n/2, 77 = 0.5, 
= kir in the case of resonance. We notice that increasing the photon number £ results in 
a decrease of the HOADS duration, which is inversely proportional to a/£I. The HOADS 
frequency is directly proportional to y^- 
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In Fig. 6, we show the time evolution of the function F±(P = 6) for 9 = n/2, r\ = 0.5 and 
(ft = kn and A = 0,1,5 in the case of one-photon transition. It is seen that the duration, 
period and strength of HOADS are much sensitive to the detuning. Especially, for a large 
detuning, for example A = 5, we observe enhanced and strong HOADS. Similar results, as 
shown in Fig. 7, are obtained for three-photon transition processes. These results show that 
large detuning can generate efficient HOADS than the resonance case. 

Finally, it would be interesting to examine the relation between the second-order ADS 
(SOADS) and HOADS. In Fig.8 and 9, we show 3D contour plots for the function F 2 (P = 2), 
9 and for 77 =| sin(2#) | /2 and r\ = 0.5, respectively, for a resonance transition. Comparing 
Fig.8 and 9 with Fig. 3 and 4, respectively, we find that there exist additional 6>-0 parameter 
regions where HOADS could occur but SOADS could not. There are also more HOADS dips 
than SOADS. For an arbitrary atomic state with r\ = 0.5, we see that HOADS with respect 
to the fluctuation of the phase is much stable than SOADS. Certainly, as demonstrated 
in our previous serial works [20-25], there exist common 9-(f> parameter regions where both 
HOADS and SOADS could be generated in the mean time. 

4. SUMMARY 

In summary, we study HOADS in a high-Q micromaser cavity by considering a two- 
leve atom interacting with an optical field through a multi-photon transition. Assuming 
that the initial atom is arbitrarily prepared and the field is initially in a vacuum state, we 
demonstrate that HOADS cannot appear if the atom is initially prepared in a chaotic state 
and that a coherent atomic state could generate less efficient and stable HOADS than an 
arbitrary one. It is found that large detuning may lead to enhanced and strong HOADS. 
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Caption of Figures 



FIG.l: The 3D contour plot of function Pi(P = 6), i] and 9 for A = 0, = kir and 

(f) = (k + 1/2)tt (k=0,l,2,...), where the dotted line is for F X {P = 6) = 0, the dotted-dashed 
lines for Pi(P = 6) = 0.01, and the dashed lines for P X (P = 6) < in the range of -0.045 to 
-0.005 with an interval of 0.005. The solid line is for r] —\ sin 9 cos 9 \ . 

FIG. 2: The 3D contour plot of function P : (P = 6), rj and 9 for A = 0, ^.r = kir and 
= kir (k=0,l,2,...), where the dotted line is for Fi(P = 6) = 0, the dotted-dashed lines for 
Fi(P = 6) > in the range of 0.01 to 0.05 with an interval of 0.01, and the dashed lines for 
Fi(P = 6) < in the range of -0.06 to -0.005 with an interval of 0.005. The solid line is for 
r] = | sin 9 cos 9 \ . 

FIG. 3: The 3D contour plot of function F 1 (P = G),<f> and 9 for A = 0, ^.t = kn and 
7] =| sin 9 cos 9 |, where the dotted line is for F\(P = 6) = 0, the dotted-dashed lines for 
Fi(P = 6) > in the range of 0.01 to 0.05 with an interval of 0.01, and the solid lines for 
Fi(P = 6) < in the range of -0.06 to -0.001 with an interval of 0.001. 

FIG. 4: The 3D contour plot of function F^P = 6), and 9 for A = 0, v^Tr = kn 
and rj = 0.5, where the dotted line is for F±(P = 6) = 0, the dotted-dashed lines for 
Fi(P = 6) > in the range of 0.01 to 0.05 with an interval of 0.01, and the solid lines for 
Fi(P = 6) < in the range of -0.06 to -0.001 with an interval of 0.001. 

FIG. 5: Time evolution of the function Pi(P = 6) for r] = 0.5, 9 = tt/2, A = and 
= kir (k=0,l,2,...): (a) f = 1 (solid line); (b) f = 4 (dotted-dashed line). 

FIG. 6: Time evolution of the function F X (P = 6) for r] = 0.5, 9 = tt/2, <p = kix 
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(k=0,l,2,...) and f = 1: (a) A = (dashed line); (b) A = 1 (dotted-dashed line); (c) A = 5 
(solid line). 

FIG. 7: Time evolution of the function Fi(P = 6) for rj = 0.5, 6 = tt/2, = kn 
(k=0,l,2,...) and f = 3: (a) A = (dashed line); (b) A = 1 (dotted-dashed line); (c) 
A = 10 (solid line). 

FIG. 8: The 3D contour plot of function F 2 (P = 2), and 9 for A = 0, y/gr = kir and 
7] =| sin 9 cos |, where the dotted line is for F 2 (P = 2) = 0, the dotted-dashed lines for 
F 2 (P = 2) > in the range of 0.01 to 0.24 with an interval of 0.01, and the solid lines for 
F 2 (P = 2) < in the range of -0.06 to -0.005 with an interval of 0.005. 

FIG. 9: The 3D contour plot of function F 2 (P = 2), and 6 for A = 0, = kit 

and rj = 0.5, where the dotted line is for F 2 (P = 2) = 0, the dotted-dashed lines for 
F 2 (P = 2) > in the range of 0.01 to 0.22 with an interval of 0.03, and the solid lines for 
F 2 (P = 2) < in the range of -0.24 to -0.01 with an interval of 0.01. 
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Xie & Rao: FIG.l/PHYSICA A 
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